Influence of Bottom Quark Jet Quenching on Single Electron Tomography of Au+Au 
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High transverse momentum single (non-photonic) electrons are shown to be sensitive to the stop- 
ping power of both bottom, 6, and charm, c, quarks in AA collisions. We apply the DGLV theory of 
radiative energy loss to predict c and b quark jet quenching and compare the FONLL and PYTHIA 
heavy flavor fragmentation and decay schemes. We show that single electrons in the pr — 5 — 10 
GeV range are dominated by the decay of b quarks rather than the more strongly quenched c quarks 
in Au+Au collisions at y/s — 200 AGeV. The smaller b quark energy loss, even for extreme opacities 
with gluon rapidity densities up to 3500, is predicted to limit the nuclear modification factor, Raa, 
of single electrons to the range Raa ~ 0.5 — 0.6, in contrast to previous predictions of Raa < 0.2—0.3 
based on taking only c quark jet fragmentation into account. 

PACS numbers: 12.38.Mh; 24.85.+p; 25.75.-q 
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Introduction. 

Recent data [1] from the Relativistic Heavy Ion Col- 
lider (RHIC) on "perfect fluidity" [2]- [5] and light quark 
and gluon jet quenching [6]- [9] provide direct evidence 
that a novel form of strongly interacting Quark Gluon 
Plasma (sQGP) is created in central Au+Au collisions 
at y/s = 200 AGeV [10]. 

In the near future, measurements of heavy quark jet 
quenching will provide further important tests of the 
transport properties of this new form of matter. In par- 
ticular, rare heavy quark jets are valuable independent 
probes of the intensity of color field fluctuations in the 
sQGP because their high mass (m c « 1.2 GeV, m& w 4.75 
GeV) changes the sensitivity of both clastic and inelas- 
tic energy loss mechanisms in a well defined way [11]- 
[17] relative to those of light quark and gluon jets [6]-[9]. 
Open heavy quark meson (D, B) tomography also has 
the unique advantage that - unlike light hadron (n, K) 
tomography that is sensitive to the large difference be- 
tween quark and gluon energy loss - gluon jet fragmen- 
tation into D and B mesons can be safely neglected. 

The "fragility" of light hadron tomography pointed out 
in Ref. [18] is primarily due to the significant reduction in 
sensitivity of the attenuation pattern to the sQGP den- 
sity when the gluon jets originating from the interior are 
too strongly quenched. In that case, the attenuation of 
light hadrons becomes sensitive to geometric fluctuations 
of the jet production points near the surface "corona" . 

Heavy quarks, especially b quarks, are predicted to be 
significantly less fragile in the DGLV [12]-[15] theory of 
radiative energy loss because their energy loss is expected 
to be considerably smaller. If radiative energy loss is 
the dominant jet quenching mechanism in the pr ~ 10 
GeV region, then heavy meson tomography could be a 
more sensitive tomographic probe of the absolute scale of 
density evolution and the opacity of the produced sQGP. 

However, one disadvantage of heavy meson tomogra- 
phy is that direct measurements of identified high pr D 
and B mesons are very difficult with current detectors 



and RHIC luminosities [19]. Therefore, the first experi- 
mental studies of heavy quark attenuation at RHIC have 
focused on the attenuation of their single (non-photonic) 
electron decay products [20]- [23]. 

Some preliminary data [24]- [25] surprisingly suggest 
that single electrons with px ~ 5 GeV may experience 
elliptic flow and suppression patterns similar to light par- 
tons. We emphasize in this letter that either result would 
have even greater implications than previously thought 
about the nature of the produced sQGP. If confirmed in 
the final analysis, the sQGP would have to be completely 
opaque to even b quark jets of pr ~ 10 GeV, in contra- 
diction to all radiative energy loss estimates so far. 

A significant complication of the heavy quark decay 
lepton measurements is that estimates in Refs. [26, 27] in- 
dicated that bottom decay leptons may in fact dominate 
electrons from charm for pr > 3 GeV in pp collisions. In 
this letter, we show that jet quenching further amplifies 
the b contribution to the lepton spectrum and strongly 
limits the nuclear modification factor of electrons in AA 
collisions. 

The preliminary electron data [24]- [25] are so surpris- 
ing that novel jet energy loss mechanisms may have to be 
postulated [28]- [31]. The elliptic flow of high pt heavy 
quarks can be accounted for, e.g., if the elastic cross sec- 
tions of all partons, including bottom, are assumed to be 
anomalously enhanced to > 20 mb, far in excess of per- 
turbative QCD predictions, up to at least pr ~ 10 GeV. 
While these enhanced cross sections could lead to heavy 
flavor elliptic flow at the pion level even at high px, they 
may greatly overestimate the attenuation of light and 
heavy flavored hadrons [31]- [33]. 

Given the critical role that single electron tomography 
of the sQGP may play in the near future, it is especially 
important to scrutinize the theoretical uncertainties and 
robustness of current predictions. This is the aim of this 
letter. 

Theoretical framework. 
The calculation of the lepton spectrum includes initial 
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heavy quark distributions from perturbative QCD, heavy 
flavor energy loss, heavy quark fragmentation into heavy 
hadrons, Hq, and Hq decays to leptons. The cross sec- 
tion is schematically written as: 
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EidMQ) 



dp 3 

D(Q — > 
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where £g) is a generic convolution. The electron decay 
spectrum, f(EtQ — > e), includes the branching ratio to 
electrons. The change in the initial heavy flavor spectra 
due to energy loss is denoted P(Ei — ► Ef). 

The initial heavy quark px distributions are computed 
at next-to-leading order with the code used in Ref. [34, 
35]. We assume the same mass and factorization scales 
as in Ref. [36] , employing the CTEQ6M parton densities 
[37] with no intrinsic kx- 

As in Ref. [15], we compute heavy flavor suppression 
with the DGLV generalization [12] of the GLV opacity 
expansion [7] to heavy quarks. We take into account 
multi-gluon fluctuations as in Ref. [8]. 

The fragmentation functions Die — > D) and D(b — > 
B), where D and B indicate a generic admixture of 
charm and bottom hadrons, are consistently extracted 
from e + e~ data [38-40]. The charm fragmentation func- 
tion [40] depends on the parameter r [41]. We take 
r = 0.04 for m c = 1.2 GcV. Bottom fragmentation in- 
stead depends on the parameter a [42] with a = 29.1 for 
nib = 4.75 GeV. The fragmentation is done by rescaling 
the quark three-momentum at a constant angle in the 
laboratory frame. 

The leptonic decays of D and B mesons are controlled 
by measured decay spectra and branching ratios. The 
spectrum for primary B — > e decays has been measured 
recently [43, 44]. The fit to this data [34] is assumed to 
be valid for all bottom hadrons. Preliminary CLEO data 
on the inclusive semi-leptonic electron spectrum from D 
decays [45] have also been fitted [34] and assumed to 
be identical for all charm hadrons. The contribution of 
leptons from secondary B decays B — > D — > e is ob- 
tained as a convolution of the D — > e spectrum with 
a parton-model prediction for b — > c decay [34]. The 
resulting electron spectrum is very soft, making it a neg- 
ligible contribution to the total, particularly at pt > 2 
GeV. The appropriate effective branching ratios are [46]: 
B(B -> e) = 10.86 ± 0.35%, B(D -> e) = 10.3 ± 1.2%, 
and B(B -> D -> e) = 9.6 ± 0.6%. 

The uncertainty in our results due to the choice of 
fragmentation and decay schemes is studied using the 
corresponding PYTHIA [47] routines, assuming Peterson 
fragmentation [48] with a range of parameters. 

To compute the medium induced gluon radiation spec- 
trum, we need to include in general three effects: 1) the 
Ter-Mikayelian or massive gluon effect [13, 14], 2) tran- 
sition radiation [49] and 3) medium-induced energy 
loss [12, 14]. In Ref. [50], it was shown that first two 
effects nearly cancel and can thus be neglected for heavy 
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FIG. 1: The differential cross section (per nucleon pair) of 
charm (upper blue) and bottom (upper red) quarks calculated 
to NLO in QCD [34] compared to single electron distribu- 
tions calculated with the fragmentation and decay scheme of 
Ref. [34]. The solid, dotted and long dashed curves show the 
effect of DGLV heavy quark quenching with initial rapidity 
densities of dN g /dy = 0, 1000, and 3500, respectively. 



quark suppression at zeroth order in opacity. We there- 
fore only compute the medium-induced gluon radiation 
spectrum [12]. We employ the effective static medium 
approximation formula 
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Here E = \Jp\ J r m 2 is the initial energy of a heavy 
quark of mass m, k is the transverse momentum of the 
radiated gluon and q is the momentum transfer to the 
jet. The opacity of the medium to radiated gluons is 
L/X g = 97ra 2 /2 / drp(r)/^ 2 (r) where n « <?(p/2) 1/3 is 
local Debye mass in a perturbative QGP. The gluon den- 
sity at proper time r is related to the initial rapidity den- 
sity of the produced gluons by p(r) ps dN g / dymR 2 with 
R — 6 fm in central collisions assuming a uniform cylin- 
der undergoing a Bjorken 1+1D expansion. Transverse 
expansion does not significantly affect the integrated en- 
ergy loss [51]. 

The Ter-Mikayelian effect at first order in opacity 
is due to a asymptotic transverse gluon mass in the 
medium, m g « n/\2. We assume a s = 0.3. The induced 
radiative energy loss fluctuation spectrum, P(Ei — > EA, 
was computed as in Ref. [8], starting from the average 
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FIG. 2: Heavy quark jet quenching before fragmentation into mesons for dN g /dy = 1000 (left) and 3500 
(right) are compared to light (u, d) quark and gluon quenching. The resulting it Raa is compared to the 
central 0-10% PHENIX data [52]. 
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FIG. 3: Single electron attenuation pattern for initial dN g /dy — 1000, left, and dN g /dy = 3500, right. The 
solid curves employ the fragmentation scheme and lepton decay parameterizations of Ref. [34] while the 
dashed curves use the Peterson function with e c = 0.06 and et = 0.006 and the decay to leptons employed 
by the PYTHIA Monte Carlo. Note that even for the extreme opacity case on the right the less quenched 
b quark jets dilute Raa so much that the modification of the combined electron yield from both c and b 
jets does not fall below ~ 0.5 — 0.6 near pr ~ 5 GeV. 



induced gluon spectrum in the effective static medium 
approximation given by Eq (1). In this approximation 
the effective static p is approximated by p((t}) with 
(t) = R/2 = 3 fm and L = R. We have checked that the 
more numerically intensive Bjorken expansion gives very 
similar results. 

Note that k max = 2x(l — x)pr in Eq. (1) instead of 
^max = xE, as in Ref. [12]. There is a 20% theoreti- 
cal uncertainty in Raa due to the range of reasonable 
kinematic bounds. 



Bottom versus Charm quark suppression 

Figure 1 shows the c and b quark distributions at 
midrapidity before fragmentation. The solid curves in- 
dicate that, at NLO, b production becomes comparable 
to c production in the vacuum only for px ~ 15 GeV. 
However, jet quenching is greater for the lighter c quark, 
and for the default gluon density, dN g /dy = 1000 [9], 
the more weakly quenched 6's dominate over the more 
strongly quenched c's for ~ 9 GeV. For more extreme 
opacities, characterized here by dN g /dy = 3500, the cross 
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over shifts down to pt ~ 7 GeV. With the fragmentation 
and decay scheme of Rcf. [34] , the electron decay distri- 
butions, c — > e and 6 — > e, are seen to cross each other at 
Pt ~ 5-5 GeV when the c and b quarks are not quenched, 
reduced to px ~ 3 GeV for dN g /dy — 3500. The electron 
results for dN g /dy — 1000, lying between the solid and 
long-dashed curves in Fig. 1, are not shown for clarity. 
Thus electrons in the pr ~ 5 GeV region are sensitive to 
b and c quark quenching. 

The parton level quenching is shown in detail in 
Fig. 2 by the nuclear modification factor, Raa(Q) = 
dNQ(p Tl dNg/dy)/dN Q (p T ,0) with Q — g,u,d,c and b. 
The left-hand side shows results for the default case, 
dNg/dy = 1000 [9] , while the right-hand side shows the 
high opacity case, dN g /dy = 3500. For comparison, we 
also show the PHENIX [52] data on the 7r° nuclear mod- 
ification factor measured in the central 0-10% of Au+Au 
collisions at = 200 AGeV. As expected, gluon quench- 
ing is largest due to its color Casimir factor and its small 
in- medium mass. The "dead cone effect" [11] is seen by 
comparing c quark quenching to light u, d quenching at 
p T < 10 GeV. For p T > 10 GeV > m c , the mass dif- 
ference between the charm and light quarks is almost 
negligible [50]. 

However, in both cases, b quark quenching remains sig- 
nificantly smaller than that of the light and charm quarks 
for pt < 20 GeV since pT/rrib is not large. The effect of 
the b mass can therefore never be neglected in the RHIC 
kinematic range. 

Figure 2 also shows an estimate of n° quenching as- 
suming 

Raa(*°) « f g Raa(9) + (1 - /„) Raa(u) , (2) 

where f g « exp[— pt/10-5 GeV] is the fraction of pions 
with a given pt that arise from gluon jet fragmentation. 
The approximate form is a fit to a leading order QCD cal- 
culation at y/s — 200 AGeV, discussed in Refs. [53, 54]. 
The approximation in Eq. (2) is strictly valid only for 
pure power law gluon and quark distributions with a 
PT-independent spectral index. However, it provides a 
simple estimate that shows that tt° quenching is primar- 
ily controlled by light quark quenching above 10 GeV. In 
addition, Fig. 2 shows that current data would be incom- 
patible with radiative g, u and d quenching if the medium 
had an opacity greater than that of the dN g /dy = 3500 
case considered on the right-hand side. 

We note that the c quark quenching predicted in Fig. 2 
with 1000 < dNg/dy < 3500 is similar to the quenching 
range predicted in Fig. 2 of Ref. [17] for the effective 
transport coefficient q = fJ- 2 /X g in the range 4 < q < 14 
GeV 2 /fm. For a c quark withp-r ~ 12 GeV, for example, 
we predict Raa{c) ~ 0.25 — 0.5 in this range, as docs 
Ref. [17] for the same factor of 3.5 variation of the sQGP 
density. 

Our primary new observation is that since b quark 
quenching is greatly reduced relative to c quenching, if 
heavy quark tomography is performed via single electron 
suppression patterns, the lower b quenching strongly lim- 



its the possible electron quenching, as we show in Fig. 3. 
For electrons arising from c fragmentation and decay, we 
again confirm the predictions of Ref. [17]. However, for 
electrons arising from b decay, there is only a modest 
amount of quenching. Note the similar magnitudes of 
heavy quark and decay electron quenching if the quark 
Pt is rescaled by a factor of <~ 2. 

In Fig. 3, the sensitivity of the electron quenching to 
variations in the heavy quark fragmentation scheme is 
shown by the difference between the solid and dashed 
curves. The solid curves are calculated as in Ref. [34] 
while the dashed curves arise when Peterson fragmenta- 
tion (e c = 0.06, eb — 0.006) is used. While there can 
be considerable differences in the fragmentation schemes 
on an absolute scale, see Fig. 4, these differences mostly 
cancel in the nuclear modification factors shown in Fig. 3. 

The yellow band corresponding to the combined c+b — > 
e electron sources shows that, in the kinematic range 
4 < pr{e) < 10 GeV accessible at RHIC, Raa(c) is domi- 
nated by b quark quenching. Even for the highest opacity, 
shown on the right-hand side, we therefore predict that 
due to the b — > e contribution 

RAA(e) > 0.5 forp T < 6 GeV . (3) 

Increasing the opacity further is not an option within the 
theory of radiative energy loss because pion quenching 
would then be over-predicted. 

The robustness of the bottom dominance in the elec- 
tron spectrum can be seen in the ratio of charm rela- 
tive to bottom decays to electrons in Fig. 4. We use the 
NLO MNR code [35] to compute heavy quark produc- 
tion for a range of mass and scale values: 1.2 < m c < 1.7 
GeV, 4.5 < m& < 5 GeV and combinations of the renor- 
malization, and factorization, fj,p, scales such that 
(Mfl/mr, H F /m T ) = (1, 1), (2,1), (1,2) and (2,2). We em- 
ploy the same (pr/tiit, hf/ttit) combinations for both 
charm and bottom to maintain the asymptotic approach 
of the distributions at high pr- In all cases, the bottom 
contribution becomes larger for pt < 5.5 GeV, even be- 
fore energy loss is applied. Changing the Peterson func- 
tion parameter, eg, from the standard values of e c = 0.06 
and £(, = 0.006 to the more delta-function like values 
of e c = €b = 10 -5 shifts the cross over to higher pt, 
more similar to the results with the FONLL fragmenta- 
tion scheme. No reasonable variations of the parameters 
controlling fragmentation of heavy quarks can make the 
bottom contribution to electrons negligible at RHIC. 

As a final check, in Fig. 5 we show the "electron 
reach" defined by the transverse momentum distribution 
of the initial heavy quarks that decay to electrons with 
Pt = 5 — 6 GeV. As can be readily seen, this range of elec- 
tron pt is sensitive to heavy quark quenching at approxi- 
mately twice this scale: pt ~ 6 — 10 GeV with hard frag- 
mentation parameters, e c = e& = 10~ 5 , and pt ^9—14 
for the standard Peterson parameters. Given the slow 
variation of heavy quark quenching in the pt ~ 10 — 20 
GeV range seen in Fig. 2, it is easy to understand why 
single electron quenching is robust to uncertainties in the 
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FIG. 4: [left] The ratio of charm to bottom decays to electrons obtained by varying the quark mass and 
scale factors. The effect of changing the Peterson function parameters from e c — 0.06, eb = 0.006 (lower 
band) to e c = = 10 _J (upper band) is also illustrated. 

FIG. 5: [right] The electron reach, defined by the distribution of initial c and b quark momenta that after 
fragmentation and decay produce an electron with pr = 5 — 6 GeV using the PYTHIA fragmentation 
scheme. 



heavy quark fragmentation scheme, as shown in Fig. 3. 
Conclusions. 

In this letter we predicted the nuclear modification fac- 
tor of single electrons, RAA(pT,mQ,dNg/dy), produced 
by fragmentation of quenched bottom as well as charm 
quark jets in central Au+Au collisions with ^/s = 200 
AGeV. We found that within the DGLV theory of radia- 
tive energy loss, b quark jets give the dominant contri- 
bution to pt ~ 5 GeV electrons, limiting Raa{g) > 0.5. 
Therefore, if the preliminary PHENIX data suggesting 
Raa{b) < 0.5 are confirmed, it will be a theoretical chal- 
lenge to devise novel energy loss mechanisms that make 



the sQGP opaque to bottom quarks of pr ~ 10 — 20 
GeV without over-predicting the observed light hadron 
quenching in the px ~ 10 GeV range. 
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